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1.  INTRODUCTION 

Numerous  systems  have  been  developed  over  the  past  25  years  for  remote  classification  of  sediments. 
The  generation  of  acoustic  impedance  as  a  function  of  depth  is  the  principle  behind  operation  of  most 
normal-incidence  acoustic  classification  systems.  A  typical  high-resolution  seismic  system  (15  kHz) 
used  to  classify  sediments  is  described  by  Lambert  et  al.  [IJ.  Profiles  of  sediment  impedance  are 
determined  from  echo  return  amplitude  and  pulse  character  using  acoustic  multilayer  theory.  Physical 
and/or  empirical  models  are  then  used  to  convert  these  in-situ  estimates  of  sediment  impedance  to 
sediment  type  or  values  of  sediment  physical  (porosity,  grain  size,  bulk  density,  permeability), 
geoacoustic  (compressional  and  shear  wave  velocity  or  attenuation,  acoustic  reflectivity)  or  engineering 
(shear  strength)  properties.  In  this  paper  we  assume  a  priori  that  a  mstic  sediment  classification 
techniques  give  accurate  estimates  of  in-situ  sediment  impedance.  We  will  instead  examine  the 
empirical  and  physical  models  that  are  used  to  estimate  sediment  type  or  values  of  sediment  properties 
from  the  in-situ  sediment  impedance. 

Perhaps  the  most  widely  used  empirical  relationships  between  sediment  acoustic  and  physical 
properties  were  developed  by  Edwin  Hamilton  in  the  1960’s  and  1970’s  [2].  Bachman  [3]  summarized 
Hamilton’s  relationships  with  an  emphasis  on  the  prediction  of  sediment  physical  properties  from 
sediment  acoustic  imp^ance.  This  approach  is  the  inverse  of  Hamiltion’s,  where  sediment  physical 
properties  are  used  to  predict  sediment  geoacoustic  properties  [4]. 

Several  factors  complicate  the  apparent  straightforward  prediction  of  values  of  sediment  properties 
from  remotely  sensed  sediment  impedance.  Firstly,  impedance  is  the  product  of  sediment  sound  speed 
and  sediment  bulk  density,  two  quantities  that  are  dependent  on  temperature,  porewater  salinity  and 
water  depth.  Because  empirical  relationships  are  based  on  measure.ments  made  at  (or  corrected  to) 
common  environmental  conditions,  it  is  imperative  that  estimates  of  sediment  properties  be  corrected 
for  differences  in  model  and  in-situ  conditions.  Secondly,  sediment  sound  speed  also  varies  with 
frequency.  These  empirical  relationships  are  usually  constructed  from  high-frequency  00-400  kHz) 
acoustic  measurements  which  can  differ  from  the  lower  frequencies  (3.5-30  kHz)  at  which  most 
remote  acoustic  measurement  systems  operate.  Depending  on  the  sediment  type,  determination  of 
frequency  dispersion  of  sediment  sound  speed  may  be  critical.  Thirdly,  both  sediment  physical  and 
acoustic  measurements  are  based  on  laboratory  measurements.  Sediment  samples  from  cores  are  often 
disturbed  during  collection,  handling  and  analysis,  resulting  in  field  measurements  that  differ  from  in- 
situ  measurements.  Fourthly,  some  empirical  relationships  are  generated  from  data  not  restricted  to 
surficial  sediments.  Impedance  vs.  grain  size  relationships  developed  from  long  core  samples  differ 
from  empirical  predictions  based  on  surficial  sediment  samples.  Finally,  variability  of  sediment 
properties  occurs  on  various  spatial  scales  and  affects  interpretion  of  acoustic  classification  data. 
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In  this  paper  new  empirical  relationships  between  sediment  impedance  and  sediment  physical 
properties  are  presented.  All  sediments  were  hand-collected  by  divers  or  carefully  sub-cored  from 
undisturbed  box  cores  to  minimize  sediment  disturbance.  The  large  number  of  sediment  samples  were 
collected  in  a  wide  range  of  environments  over  the  last  12  years.  All  data  are  based  on  samples 
collected  from  the  upper  0.5  m  of  sediment.  The  in-situ  effects  of  sediment  temperature,  salinity  and 
water  depth  are  eliminated  by  the  introduction  of  an  environmentally  independent  index  of  impedance 
(lOI).  The  effects  of  frequency-dependent  impedance  will  be  the  subject  of  future  papers. 


2.  EXPERIMENTAL  AND  ANALYTICAL  METHODS 

Geoacoustic  and  physical  property  measurements  were ;  made  from  sediments  collected  in 
polycarbonate  plastic  cores  that  were  6. 1-cm  in  diameter  and  45-cm  long.  Each  core  was  bevelled  at 
one  end  to  facilitate  the  careful  manual  penetration  into  the  sediment.  Cores  were  capped  at  both  ends 
immediately  upon  collection  to  retain  the  water  overlying  the  sediment  and  kept  in  an  upright  position 
during  transport  to  the  laboratory  for  analysis.  Collection,  measurement,  and  handling  procedures 
were  designed  to  minimize  sampling  disturbance  and  to  maintain  an  intaci  sediment- water  interface 
within  the  core  samples.  The  cores  penetrated  to  a  depth  of  9  to  40  cm  into  the  sediment,  depending 
on  the  sediment  texture.  Of  the  11  experimental  sites,  four  sites  were  sampled  from  boxcores 
(Montauk  Point,  Quinault  Range,  Arafura  Sea,  and  Russian  River)  and  the  rest  were  sampled  directly 
by  divers  [5-13]. 

Measurement  of  sediment  sound  speed  was  made  within  24  hours  of  collection,  once  the  sample&  had 
equilibrated  with  laboratory  temperature.  Sediment  sound  speed  and  attenuation  were  measured  at  1- 
cm  intervals  using  a  pulse  technique  [8].  Time  delay  measurements  of  a  400-kHz  cw  pulse  were  made 
on  cores  through  the  sediments  and  a  distilled  water  reference  using  an  Underwater  Systems  model 
USI-103  transducer-receiver  head.  Differences  in  time  delay  between  distilled  water  and  sediment 
cores  were  used  to  calculate  sediment  sound  speed.  Sediment  compressional  wave  attenuation  was 
calculated  as  20  log  of  the  ratio  of  received  voltage  through  distilled  water  to  received  voltage  through 
sediment.  Sound  speeds  were  corrected  to  a  common  temperature,  salinity  and  pressure  (23'C,  35%o, 

I  atm)  after  Hamilton  [14]. 

Samples  were  extruded  from  sediment  cores  upon  completion  of  acoustic  measurements  and  sectioned 
at  2-cm  intervals  (1-cm  intervals  in  cores  from  L.I.  Sound)  to  determine  sediment  porosity  and  grain 
size  distribution.  Porosity  was  determined  from  weight  loss  of  sediment  dried  at  105'C  for  24  hours. 
Sediment  grain  size  was  determined  from  disaggregated  samples  by  dry  sieving  for  sand-sized  particles 
and  by  use  of  a  Micromeritics  Sedigraph  for  silt-  and  clay-sized  particles  when  samples  were  collected 
from  muddy  environments  [10]. 

Values  of  sediment  sound  speed  are  expressed  as  the  ratio  of  measured  sediment  speed  to  measured 
speed  of  the  overlying  water  in  the  core  (same  temperature,  salinity,  and  depth).  Values  of  attenuation 
are  expressed  in  units  of  dB  m  'kHz  ’  and  are  identical  to  the  constant,  k,  described  by  Hamilton  [15]. 
Sediment  impedance,  Z,  is  calculated  as  the  product  of  the  measured  values  of  sediment  sound  speed 
and  the  calculated  values  of  sediment  density  and  has  the  units  of  g  cm'^  s  ‘.  Density  is  calculated  from 
sediment  porosity,  grain  density  and  seawater  density  [8].  Values  of  the  index  of  impedance  (lOI)  are 
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calculated  as  the  product  of  the  sediment  velocity  ratio  and  the  sediment  density.  The  ratio  is 
unitless,  therefore  the  lOI  has  units  of  density  (g  cm  *). 


3.  RESULTS  AND  DISCUSSION 

Sediment  physical  and  geoacoustic  properties  were  measured  on  211  cores  collected  at  11  different 
sites  worldwide.  Table  1  lists  the  locations  from  which  surficial  sediment  cores  were  collected,  the 
sediment  characteristics  including  the  average  measured  values  of  sound  speed  (V^  ratio),  sound 
attenuation  (k),  porosity  and  mean  grain  size  {4>)  and  the  average  calculated  values  of  impedance  (Z) 
and  index  of  impedance  (lOI)i  The  number  of  data  pairs  ratio  and  porosity)  acquired  in  order  to 
calculate  impedance  values  is  tabulated  as  the  number  of  data  points  at  each  location.  This  number 
represents  the  frequency  with  which  porosity  values  were  measured  at  the  same  core  depth  interval 
as  sediment  sound  speed  values. 

As  suggested  in  an  examination  of  similar  data  [13],  sediment  physical  and  geoacoustic  properties  are 
quite  variable  on  scales  of  a  kilometer  or  less.  The  large-scale  variability  of  sediment  properties  is 
exhibited  by  the  Long  Island  Sound,  Mission  Bay,  St.  Andrew  Bay  and  La  Spezia  locations.  These 
locations  are  divided  into  sedimentary  provinces  where  a  particular  sediment  type  predominates. 
Small-scale  variability  of  sediment  properties  within  a  particular  sediment  type  is  demonstrated  in 
Table  2.  The  coefficient  of  variation  of  the  ratio,  attenuation  (k),  porosity  mean  grain  size  ((p),  bulk 
density  (p)  and  lOI  is  calculated  as  the  standard  deviation  divided  by  the  mean  and  expressed  as  a 
percentage.  The  coefficient  of  variation  (CV)  of  impedance  (Z)  is  identical  to  the  CV  of  the  101 
because  the  CV’s  of  the  components  are  the  same.  Several  general  observations  can  be  made 
pertaining  to  this  variability.  Sound  speed  (or  velocity  ratio)  is  the  least  variable  sediment  property; 
compressional  wave  attenuation  is  the  most  variable.  However,  in  locations  such  as  Jacksonville  (II), 
St.  Andrew  Bay  and  La  Spezia  (Tellaro)  mean  grain  size  exhibits  great  variability  due  to  occurrence 
of  discrete  layers  of  different  sediment  types.  The  variability  of  sediment  impedance  (or  lOI)  appears 
correlated  with  the  variability  of  sediment  bulk  density  rather  than  sound  speed.  Overall  variability 
of  all  sediment  properties  is  lower  at  muddy  rather  than  sandy  sites.  The  greatest  variations  are  found 
at  sandy  sites  that  contained  considerable  amounts  of  shell  material  or  at  sites  that  contained  mixtures 
of  sand  and  mud. 

The  variability  of  sediment  physical  and  geoacoustic  properties  presented  in  Table  2  was  determined 
from  replicate  sediment  samples  collected  from  sites  generally  within  an  area  with  a  100-m  radius. 
Except  for  the  La  Spezia  and  Panama  City  (II)  data,  replicate  cores  were  collected  within  a  short 
(usually  a  week)  period  of  time.  Coefficients  of  variation  calculated  for  cores  collected  within  a  week 
of  each  other  at  the  La  Spezia  and  Panama  City  (II)  sites  were  similar  to  the  longer  term  data  sets. 
We  therefore  suggest  that  the  CV’s  presented  in  Table  2  represent  variability  that  can  be  expected  at 
scales  used  to  characterize  sediment  physical  properties  remotely  by  acoustic  methods.  Given  the  high 
variability  of  compressional  wave  attenuation,  precise  prediction  of  this  geoacoustic  property  may  be 
difficult. 
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Table  1 .  Data  sources  and  mean  values  of  geoacoustic  and  physical  properties. _ 

Experimental  Site  Water  Location  Sediment  No.  of  Vp  Atten.  Porosity  (p  Z  lOI 

Depth  Type  Data  Ratio  (k)  (%)  (gem’s)  (gem’) 

(m)  Pts.  (X  ItP) 
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Table  2.  Coefficients  of  variation  (%)  for  sediment  physical  and  geoacoustic  properties. 


Experimental  Site 

Vp 

Ratio 

Atten. 

m 

Porosity 

P 

lOl 

Long  Island  Sound 

NWC 

0.4 

— 

1.6 

2.4 

1.5 

1.4 

FOAM 

0.8 

— 

7.5 

11.8 

6.3 

6.8 

Mission  Bay,  CA 

fine 

1.1 

14.8 

— 

10.8 

— 

— 

coarse 

1.0 

25.7 

— 

7.4 

— 

— 

Montauk  Point,  NY 

1.0 

16.6 

3.2 

3.7 

1.0 

1.4 

Quinault  Range,  WA 

1.2 

43.6 

5.1 

3.5 

1.8 

3.1 

Charleston,  SC 

1.0 

37.9 

4.8 

19.1 

1.6 

1.9 

Arafura  Sea,  AUSTRALIA 

C.4 

34.7 

5.8 

14.9 

4.5 

4.6 

Panama  City,  FL 

I 

0.9 

16.4 

3.1 

4.2 

1.0 

2.0 

n 

1.6 

30.5 

6.3 

12.1 

2.7 

4.1 

Jacksonville,  FL 

I 

1.0 

31.3 

4.1 

6.1 

1.3 

2.1 

n 

1.8 

n.i 

7.8 

81.6 

2.6 

3.3 

St.  Andrew  Bay,  FL 

f.sand/clay 

4.9 

71.6 

28.4 

58.3 

19.8 

24.2 

clay 

0.1 

24.9 

2.0 

3.7 

2.3 

2.3 

fine  sand 

0.3 

25.5 

2.1 

3.5 

0.6 

0.8 

Russian  River,  CA 

0.5 

14.6 

6.2 

7.1 

3.8 

4.3 

La  Spezia,  ITALY 

Santa  Teresa 

1.5 

80.3 

13.0 

9.8 

9.2 

11.5 

Portovenere 

0.4 

35.4 

6.2 

2.1 

4.4 

4.3 

Venere  Azzura 

1.9 

23.8 

5.6 

35.2 

2.2 

4.3 

Tellaro 

4.2 

39.8 

9.8 

36.6 

4.3 

8.0 

Monasteroli 

2.6 

15.3 

6.8 

27.6 

1.9 

3.7 

Diga 

0.2 

88.9 

2.4 

0.8 

1.8 

1.8 

Viareggio 

0.6 

51.2 

7.5 

2.2 

4.6 

4.9 

Polynomial  or  linear  regressions  are  constructed  for  the  various  geoacoustic  and  physical  property 
measurements  made  at  all  locations  in  Table  1.  The  regression  equations,  coefficient  of  determination 
(r^)  and  the  F-value  for  the  regression  equations  are  tabulated  in  Table  3.  The  coefficient  of 
determination  is  used  to  indicate  the  proportion  of  variation  of  one  variable  determined  by  the 
variation  of  the  other.  In  every  case  in  Table  3  the  calculated  F-value  shown  in  Table  3  exceeds  the 
tabulated  F-value  with  appropriate  degrees  of  freedom  by  more  than  a  factor  of  four,  indicating  that 
the  regression  generates  a  "truly  satisfactory  prediction"  [3J. 
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FIG  2.  Relationships  between  geoacoustic  and  physical  properties  and  lOI. 

Data  are  concentrated  in  two  regions  in  the  regression  plots  corresponding  to  sands  and  muds.  The 
intervening  data  are  widely  dispersed  and  correspond  to  various  mixtures  of  coarse  and  fine  sediments 
(Figs.  1,  2).  Data  dispersion  due  to  over-consolidated  clays  are  indicated  in  the  plots  by  fine-grained 
sediments  (high  4>  values)  with  low  porosity  or  high  bulk  density  values.  Values  of  ratio  are 
regressed  against  values  of  physical  and  geoacoustic  properties  and  displayed  in  Table  3,  but 
regressions  predicting  sound  speed  in  m  s'*  are  obtained  by  multiplying  the  regression  coefficients  by 
the  the  sound  speed  in  water  for  the  particular  environment  of  interest.  A  curious  similarity  in  the 
relationships  between  the  vs.  bulk  density  data  and  the  vs.  lOI  data  is  noticeable  in  Figs.  1  and 
2.  The  similarity  is  apparent  in  the  mean  grain  size  v^.  bulk  density  data  and  the  mean  grain  size  vs. 
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lOI  data.  This  is  a  result  of  the  fact  that  the  variations  in  bulk  density  values  exert  a  strong  influence 
on  calculated  values  of  lOI.  Attenuation  (k)  is  not  correlated  with  mean  grain  sixe  (Fig.  1)  or  any 
other  geoacoustic  and  physical  property.  Empirical  relationships  between  sediment  physical  and 
geoacoustic  properties  exist,  but  these  relationships  are  only  correlative,  not  causal. 


Table  3.  Regression  equations  and  statistics. 


Regression 

Variables 

Regression 

Equation 

r" 

F  value 

Vp  Ratio  vj.  ^ 

Vp  - 

=  1.180  -  0.034<^  +  0.0013^ 

0.820 

2378.27 

Vp  Ratio  vj.  n  (%) 

Vp  - 

=  1.574  -  0.015/1  +  O.OOOl/i’ 

0.954 

11710.7 

Vp  Ratio  vs.  p  (g  cm  ’) 

Vp- 

=  1.623  -  0.936p  +  0.3417p’ 

0.944 

9629.79 

Vp  Ratio  VJ.  Z  X  10’  (g  cm’ s) 

Vp  - 

=  1.174  -0.2072-0.05602? 

0.972 

19634.0 

4>  vs.  lOl 

4  = 

20.23  -  9.48/0/  +  0  667/0/’ 

0.828 

2116.18 

n  (»)  vs.  lOl 

n  = 

202.14  -  120.70/0/  +  21.598/0/’ 

0.996 

142632 

p  (g  cm  ’)  vs.  lOl 

P  = 

-0.502  +  1.802/0/ -0.3050/0/’ 

0.996 

126646 

Vp  Ratio  vs.  lOl 

Vp- 

=  1.173  -  0.315/0/  +  0.1296/0/’ 

0.972 

19708.0 

vs.  n  (%) 

-4.55  +  0.169/1 

0.805 

3792.13 

<t>  vs.  p  (g  cm  ’) 

22.85  -  10.275p 

0.809 

3920.22 

The  Index  of  Impedance  is  a  \  cry  good  predictor  of  sediment  sound  speed,  bulk  density,  porosity  and 
mean  grain  size.  The  values  for  the  regressions  in  Fig.  2  range  from  0.805  to  0.996.  Hence,  lOI 
can  be  used  with  confidence  to  predict  sediment  type.  The  over-consolidated  clays  separate  from  the 
main  trend  of  the  data  in  the  mean  grain  size  vs.  lOI  relationship  as  a  field  of  high  high  lOI  values 
(Fig.  2).  The  apparent  uniqueness  of  these  data  is  significant  to  the  prediction  of  sediment 
characteristics  from  impedance  measurements. 

Sediment  compressional  wave  attenuation  (k)  is  not  predicted  satisfactorily  by  lOI  (the  regression  is 
not  shown  in  Fig.  2).  Some  or  all  of  the  following  explanations  may  apply  for  a  lack  of  correlation 
between  attenuation  and  lOI.  Firstly,  attenuation  is  a  highly  variable  parameter  of  shallow-water 
sediments  and  is  therefore  notoriously  difficult  to  predict.  Secondly,  the  proportion  of  attenuation  due 
to  scattering  of  acoustic  energy  is  large  in  relation  to  the  attenuation  due  to  friction  (i.e.,  intrinsic 
attenuation)  when  the  grain  sizes  are  large  compared  with  the  acoustic  wavelength.  This  effect  causes 
a  wide  dispersion  of  measured  attenuation  values,  especially  when  shells  are  a  constituent  of 
sediments)  [9,10,13].  Thirdly,  attenuation  is  strongly  affected  by  any  sediment  disturbance  which  may 
occur  during  sampling  or  handling.  Finally,  the  attenuation  measurement  techniques  on  small  (6. 1-cm 
diameter)  cores  may  be  called  into  question. 
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Relationships  between  impedance  and  various  sediment  parameters  are  similar  in  trend  and  dispersion 
to  the  relationships  depicted  in  Fig.  2  for  lOI.  The  impedance  values  in  Table  1  and  the  equation  in 
Table  3,  however,  were  calculated  from  data  that  were  adjusted  to  a  single  temperature,  salinity  and 
depth.  The  adjustment  to  the  sound  speed  data  accounts  for  the  similarity  in  the  values  for 
impedance  and  lOI  regressions.  The  significance  in  using  lOI  instead  of  impedance  for  predictions  is 
accentuated  in  environments  where  there  is  a  large  seasonal  change  in  bottom  water  conditions. 
Different  values  in  impedance  from  a  particular  area  result  from  seasonal  variation  in  temperature  and 
salinity.  For  example,  a  fine  sand  from  La  Spezia,  Italy  impedance  varies  from  3.22  x  10*  g  cm  ^s  ‘ 
in  winter  to  3.27  x  10*  g  cm'^s  *  in  summer.  A  mud  (silty  clay)  from  Long  Island  Sound,  USA  varies 
from  2.00  x  10*  g  cm'^s  *  in  winter  to  2.06  x  10*  g  cm'^s  '  in  summer.  By  contrast,  values  of  lOI  for 
surface  sediments  from  La  Spezia  and  Long  Island  Sound  are  2. 145  and  1 .376  g  cm  *,  respectively, 
regardless  of  season. 

The  de  facto  correction  of  impedance  values  for  environmental  conditions  using  sound  speed  ratio 
values  instead  of  in-situ  sound  speed  values  does  not  significantly  improve  predictions  for  different 
sediment  types.  In  comparison  to  naturally  occurring  variability  of  sediment  properties,  the  magnitude 
of  the  correction  is  sm^l.  In  fact,  the  same  amount  of  data  dispersion  around  the  predicted  regression 
line  occurs  in  plots  of  impedance  values  and  lOI  values. 

Empirical  relationships  among  sediment  physical  and  geoacoustic  properties  presented  here  are  similar 
to  those  reported  by  Bachman  [3]  for  continental  terrace  sediments.  Predicted  sound  speed  in  sands, 
however,  were  lower  (25  to  75  m  s  ‘)  than  those  predicted  by  Bachman’s  empirical  relationships.  Both 
Bachman’s  and  our  empirical  relationships  yield  approximately  the  same  sound  speed  predictions  for 
muddy  sediments.  As  expected  for  sandy  sediments,  the  acoustic  impedance  (Z)  estimates  reported 
here  are  slightly  lower  than  predictions  made  by  Bachman.  Therefore,  for  a  given  measured  acoustic 
value  of  impedance  or  lOI  our  relationships  indicate  that  sandy  sediments  have  slightly  lower  porosity, 
higher  bulk  density  and  smaller  mean  phi  size  (larger  grain  diameter)  values  than  predicted  from 
Bachman’s  relationships. 
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